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A common philosophy in control theory is the control of disorder by order. Control of decoherence is no 
exception; strategies aimed at suppressing quantum decoherence adopt this point of view. Here we predict 
an anomalous phenomenon in open quantum systems-control of disorder by (even more) disorder. It is 
shown that suppression of decoherence can be achieved using the most disordered white noise field, 
specifically a white Poissonian noise field. This phenomenon seems to be another anomaly in quantum 
mechanics and may offer a new strategy in quantum control practices. 

Decoherence is the deterioration of quantum information in a system due to inevitable interactions with the 
environment or bath 1 " 5 . Suppression of decoherence is one of the paramount challenges in quantum 
control practices and requires accurate control of the system dynamics. Here we predict an anomaly: 
suppression of decoherence can be be achieved using uncontrollable white noise fields. By increasing the strength 
of noise signals, a two-level system becomes less coupled to its environment and even remains in the initial state 
for a period of time. The aberrant effect reveals a new physical mechanism in quantum control theory, and in 
practice may offer the possibility of control by uncontrollable white noise. We term the phenomenon as control of 
decoherence with no control. 

Noise is a source of disorder. White noise, whose spectrum has equal power within any equal interval of 
frequencies, is the extreme of disorder (in comparison with coloured noise). Over a decade ago, people began to 
notice in classical systems that noise leads not only to nuisance but also to advantages. A remarkable example is 
that an external coloured noise can suppress the intrinsic white noise 6 ' 7 . While it is a surprise, the phenomenon fits 
well with the common philosophy - control of disorder by order. This philosophy has been carried out in classical 
noise control and extended to suppression of decoherence in quantum dynamical processes. External field control 
of quantum decoherence dates back to the spin echo technique 8 . This technique is used to suppress the inhomo- 
geneous spin dephasing by applying a n inversion pulse and has been developed to tackle general decoherence 9 " 17 . 
However, the philosophy, control of disorder by order, remains the same for both the classical and the quantum 
mechanical processes. Now the predicted anomaly is opposite to the common philosophy; the most disordered 
white noise is used to control less disordered decoherence, which is characterized by a quantum stochastic process 
with coloured noise or by a non-zero correlation function over finite time. Seeing that the setting in use is 
exclusively quantum mechanical, it appears that the phenomenon is another anomaly in quantum systems. 

Quantum mechanically, the dynamical process of a system plus its environment is governed by the total 
Hamiltonian, 

H tot = Hs(f)+H B + HsB, (1) 

where H s (t) and H B are the system Hamiltonian, embedded with white noise, and the environment Hamiltonian, 
respectively. The system dynamics is normally characterized by master equations. The system-bath interaction 
H SB is the source of decoherence. 

Results 

Protocol of decoherence suppressions with no control. We now introduce our protocol of decoherence 
suppressions, in particular suppression of dissipation. Dissipation is a decoherence process caused by the 
exchange of energy between a system and the environment. Provided that a two -level system is in its 
excitation state, the environment induces the system to give off energy and to decohere to the ground state 
[Fig. 1(a)]. If nature happens to have the white noise which we have required above, the decoherence can be 
suppressed spontaneously [Fig. 1 (b)] . This required noise is described by white noise c(t) = n(J, W,t), in particular 
the biased Poissonian white noise with the strength J 18 . We name the average time interval between two neighbour 
noise signals as l/W = Tin, where T is a time scale, different for variant systems, and n is the noise arrival 
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number 18 . When II W goes to zero, c(t) corresponds to the con- 
tinuous-time white noise process, where / is the only parameter. If 
1/W is finite, c(t) can be the biased Poissonian white shot noise, 
which is essentially different from well -controlled, idealized or 
non-idealized, pulse sequences. We can tune the parameter W 
towards the continuous time limit. In what follows, we will study 
the system responses to the noise c(t). 

Fidelity preserved by white noise. Consider a dissipative model for 
the two-level system, described by a non-Hermitian Hamiltonian in 
the exact quantum Stochastic Schrodinger equation [See Method] 19 ' 20 , 

H ss ( t) = [co + c( t) } o z 1 2 + iz*ga _ - igQo + a _ , (2 ) 

where co is the bare-energy spacing and g is the coupling strength 
between system and environment. c(t) is the above-mentioned white 
noise signal and Q satisfies a nonlinear differential equation 
Q(t)=gy/2+[ — y + ico + ic(t)]Q+gQ 2 , with a boundary condition 
Q(0) = 0 17 . The correction function of this process is G(t 9 s) = 

Z e -y\ f - s \ m Here y characterizes the environmental memory in the 

Ornstein-Uhlenbek process and is inversely proportional to the 
environmental memory time. The values of y can be used to 
somehow determine the degree of non-Markovianity. The larger y 
is, the more Markovian the environment is. y —> °o corresponds to 
the white noise model and indicates the Markov limit. This colored 
noise is formulated as z* +At = z* — yz*At + ^/At/2yw*, where w* is a 
complex Wiener process 17 . 

Suppos e that the system initial state is \i// 0 ) = The fidelity 
T(t) = \J ('AoIpJ'Ao)' qualifying the survival probability, evolves 
according to 

T{t) = e~& m(m , (3) 

where is the real part of the input function 17 . Below we will 
numerically study the noise effects on fidelity during time courses. 

Figures 2 show F(t) vs time for y = 0.2 and 0.5, subject to different 
c(t). Suppression of dissipation is excellent in all cases with bigger W, 
in particular for less Markovian environments. Different values of W 



represent different physics. Smaller W corresponds to a sequence of 
noise shots with random amplitudes and sparser random arrival 
moments, as illustrated by W = 200/ T. It suppresses dissipation to 
some extent, but not as efficient as bigger W. When W is big enough, 
all different values of W tend to result in the same fidelity and give an 
identical curve, as exemplified by the red solid curves with W = 
1000/T. These curves correspond to the continuous-time white 
noise. 

The quality of suppression depends mostly on the strength / of the 
continuous -time white noise and environmental non-Markovianity 
y. The two figures show clearly that the larger / is, the better the 
quality of suppression is. The parameter y seriously influences the 
quality of suppression as well. Suppression becomes worse at y = 0.5 
than at y = 0.2, and even worse when the environment is more 
Markovian. Eventually, suppression is invalid in a complete Mark- 
ovian environment (y — » °o), corresponding to quantum white noise. 
Significantly, it shows that white noise cannot suppress white noise, 
i.e. complete disorder cannot suppress complete disorder. 

Now we look into the detailed roles that the parameters / and W 
play. Figure 3 plots the fidelity contour as a function of /and Wat two 
time moments t = 50 Tand t = 100 T. The regions where T> 0.99 
are highlighted. The fidelity is saturated at W > 600/ T for both 
figures, where a discrete random pulse sequence becomes the con- 
tinuous-time white noise. While the fidelity is excellent for bigger 
values of /, it seems not to be saturated with /. The bigger / is, the 
better the fidelity is. 

The numerical results presented by figures could be valid for vari- 
ous physical systems with corresponding characteristic values of T. 
For example, co is approximately 10 9 -10 10 Hz in a superconducting 
flux qubit 21 . The relaxation time is T x ~ 1 ^s such that the time scale 
T ~ 5 ns and the dimensionless coT ~ 5 as taken in Fig. 2. The 
required noise strength / should be more than 10 9 Hz in order to 
successfully suppress decoherence. 

Discussion 

The perfect suppression could be justified by the following argument. 
By integrating over equation (5) (see Methods), for a long time limit, 
one can always write, 
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(a) 7^0.2 (b) y=0.5 




t/T t/T 

Figure 2 | ) vs time. The environmental memory parameter y is taken as 0.2 and 0.5 respectively. We compare the free dynamics (green-solid curves) 
with those subject to different noises, red curves to W= 1 000/ T and blue to W= 200/ T. The solid, dashed and dotted curves represent different strengths 
/ = 15co, So and / = 3oj respectively. Here coT = 5 and g = OAco. 



(*Aol<Aa 



-Si 



Nap(t)hap(t)dt 9 



(4) 



where a's denote a set of complete bases. Here N(t) and /z(f) represent 
a noisy matrix and a system dynamical matrix (see example in 
Methods). If each nonzero element of AT is a fast- varying noise func- 
tion of time and the corresponding element in h is much slower and 
weaker, the integral of each term N a p(f)h x p(t) could be zero. In our 

model, N(t)cce~ l Jo c<yS ^ ds and h(t) = —Q(t){\j/ 0 \\j/ t ) are c-numbers. 
For bigger values of J,h(t) is slower and weaker than N(t). When h(t) 
is such slow in comparison with N(t) that it can be treated as time 
independent, it is easily to prove J 0 °° N(t)h(t)dt = 0 by using the 
properties of a continuous-time white noise. 

The signal c(t) mimics natural white noise and is also associated 
with dephasing processes. Our results demonstrate that the existence 



of noise c(t) (or dephasing) significantly inhibits the dissipation. This 
well explains the reality where the relaxation time T x is longer than 
the dephasing time T 2 for all systems; dephasing has suppressed 
dissipation spontaneously. 

It is easy to discriminate our approach from the well- discussed 
dynamical decoupling, e.g. 11 " 13 , since the latter control method is 
realized through designed sequences of pulses in frequency or the 
arrival time whereas white noise is ubiquitous and beyond artificial 
control. Yet it is also dramatically different from the stochastic res- 
onance 22 " 24 , where a dynamical system is subject to both periodic 
forcing (for classical systems) or Floquet Hamiltonian (for quantum 
systems) and random noise may show a resonance or coherent beha- 
vior which is absent when either the periodical driving or the per- 
turbation is absent. The stochastic resonance is achieved when the 
noise amplitude is optimized so that it is not too weak to attain the 
threshold of the system signal detector and is not too strong to 



(a) t=50T 



(b) t=100T 
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Figure 3 | Parameter thresholds and effective regions at two time moments. Here coT = 5 and g = OAco. Fidelity where T > 0.99 is distinguished. The 
non-Markovian environment is characterized by y = 0.2. 
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overwhelm the system coherent behavior. On the contrary, our con- 
trol works when the amplitude of the white noise / is above a lower 
bound. The stronger noise is, the better the control is. More precisely, 
all known control approaches fit well with the philosophy-control of 
disorder by order, while ours, control of disorder by even more 
disorder, is opposite. 

Methods 

Consider the system-bath interaction H SB = LB ] + VB for simplicity. The exact 
stochastic Schrodinger equation is 19 - 20,25 : 



id t ilf t (z*)=H ss (t)il/ t (z*)- 



(5) 



where H ss = H$(t) + iLz* — iL^O is an exact system Hamiltonian and L = go- for our 
two-level system. O is a combination of system operators and environmental noises 
satisfying consistency conditions 19 . Each quantum trajectory is accompanied by a 
special process z*, and the system density matrix is given by p t = M[ t (z*))(il/ t (z*) \ ] . 

Suppose that the system is initially at |i^ 0 > = a vector in the completed set | v)'s. 
The system Hamiltonian and the coupling operator could be generally expressed by 
H s (t) = E v tt> v (£)|v)(v| and! = S^ v C^ v |/^><v| respectively, 



h a p = z*Cp a 3 Pll ((x.\\l/ t )- ^C^J^6^(V\0\^)^i'\\l/ t ) 



(6) 



and No,r = e 1 X m ^] d \ The white noise is embedded in the differences co ff 



wp.T = \J (i/f 0 |M[|i/^)(i/^|]|i/f 0 ) is obtained by an ensemble average over the integral 
results in Eq. (4). For the two-level system initially at 1 1), we can specifically write Eq. 
(4) as, 



(^ol^oo>-l = - f Q(t)N*(M^t) 
Jo 



dt 



(7) 



whereN(0 = e- ! i"o [£U+c(s)] " 5 . 

The paper employs the biased Poissonian white noise 18,26 , c(t) = HjXj3(t — tj) 
satisfying the following statistical properties 

M[c(t)]=JWJ=M[xj], 

M[c(t)c(s)} -M[c(t)]M[c(s)} = WM \xf\ S(t-s), 



(8) 



where xjs are noise heights. Details of numerical realization of the noise can be found 
in Ref. 26. 
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